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We have studied thermal gradients in thin Cu and AIMn wires, both exper- 
imentally and theoretically. In the experiments, the wires were Joule heated 
non-uniformly at sub-Kelvin temperatures, and the resulting temperature gra- 
dients were measured using normal metal-insulator-superconducting tunnel 
junctions. The data clearly shows that even in reasonably well conducting 
thin wires with a short f~ 10 fim) non-heated portion, significant tempera- 
ture differences can form. In most cases, the measurements agree well with 
a model which includes electron-phonon interaction and electronic thermal 
conductivity by the Wiedemann- Franz law. 

PACS numbers: 73.23.-b, 72.10.Di, 74.50.+r 



1. INTRODUCTION 

Temperature is naturally the most critical quantity in studies of any 
thermal properties of materials, and its measurement is typically a non- 
trivial task. In nano- and mesoscopic structures at low temperatures, this 
task is made even harder by the small size of the samples and their sen- 
sitivity to any external noise power, which is caused by weakness of the 
electron-phonon interaction.^ Because of this weakness, electrons can be eas- 
ily overheated (or cooled) with respect to phonons at temperatures below 
IK, a phenomenon known as the hot-electron effectP In this quasiequilib- 
rium regime, valid in many situations, the electronic and phononic degrees 
of freedom attain their own internal equilibrium temperatures, with a power 
flowing between them. Given this regime, it is possible to measure the elec- 
tron temperature, if a suitable thermometer is found. 
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A good choice for electron thermometry below IK is normal metal- 
insulator-superconductor (NIS) tunnel junction thermometerspE' known for 
their sensitivity. They can also be fabricated to a very small size ~ 100 x 100 
nm, so that it is possible to measure local temperatures. However, there are 
few studies describing how sample geometry affects temperature profiles in 
mesoscopic samples, and consequently how it complicates the interpretation 
of temperature measurements. 

In this paper, we discuss experimental and theoretical results on temper- 
ature gradients in mesoscopic metal wires with non-uniform Joule heating. 
The experiments were performed with NIS junctions, therefore allowing us 
to measure temperatures at several locations. We observe that there are 
significant temperature gradients in copper and aluminum manganese wires, 
even if the non-heated portion of the wire is only ~ 10 ^m long. The exper- 
imental observations can be explained successfully by numerical solution of 
a non-linear differential equation describing the heat balance, which incor- 
porates the Wiedemann-Franz law and electron-phonon coupling theory. In 
addition, we also present numerical results on temperature profiles caused 
by non-uniform cooling, modeling e.g. tunnel junction coolers.^ 

The paper is organized as follows: Sec. 2 discusses the theory. In Subs. 
2.1 we briefly review the theory of electron-phonon interaction in disordered 
films, followed by the numerical results on non-uniform Joule heating in 2.2 
and on cooling in 2.3. Section 3 presents the experimental techniques, with 
the experimental results presented in Sec. 4. Conclusions are drawn in Sec. 
5. - 



2. THEORY AND NUMERICAL RESULTS 

2.1. Electron-phonon coupling in metals 

Electron-phonon (e-p) scattering is a critical process for understanding 
how electron gas relaxes energy. It is the dominant mechanism for electron 
energy loss below 1 K, as photon radiative losses are very small except in 
the smallest samples!^ The strength of e-p coupling depends significantly 
on several factors: the material in question; the level of disorder in the metal, 
parametrized by the electron mean free path Zj^l an d the type of scattering 
potentialP 

In general, the electron-phonon scattering rate has a form 

— = aT?, (1) 
and the corresponding net power density transferred from hot electrons to 



Influence of Temperature Gradients on Tunnel Junction Thermometry below 1 K 



phonons is 

p = X( T n_ T n )i (2) 

where £ describes the strength of electron-phonon coupling, T e is the elec- 
tron and T p the phonon temperature, and n = m + 2. The exact form 
of coupling constant £ and the exponent n is determined by the disorder, 
mainly depending on the parameter ql, where q is the wavevector of the 
dominant thermal phonons. 

In ordered metals, defined as ql > 1, electrons can scatter either from 
longitudinal only, or from longitudinal and transverse phonons depending 
on temperature and material.^ When scattering happens only from lon- 
gitudinal phonons, the temperature depende nce f or scattering rate r~_ p is 
m = 3 and, for the heat flow in Eq. (HJ, n = 5JIHU In this case, the coupling 
constant £ is only a material dependent parameter. If electrons interact 
dominantly with transverse phonons, m = 2, n = 4 and £ cx 1/Z. 9 However, 
if the scattering rates for transverse and longitudinal phonons are are of 
the same magnitude, m can vary continuously between m = 2 — 3, so that 
n = 4 - 5 and S oc l' 1 - 1°. m 

In disordered metals, where ql < 1, electrons scatter strongly from 
impurities, defects and boundaries, and the situation is more complicated to 
model because of the interference processes between pure electron-phonon 
and electron- impurity scattering events. However, there is a theory that 
includes electron- impurity scattering by vibrating and static disorder Pin the 
case of fully vibrating impurities (following the phonon mode) m = 4, n = 6 
and X cx I. On the other hand, if the scattering potential is completely static, 
electrons interact only with longitudinal phonons and at low temperatures 
m = 2, n = 4 and S oc l/l. Between these two extremes there is a region, 
where the scattering potential is a mixture of the two, and theory predicts 
exponents m ranging between 2-4 and S depending on I as l k , where k ranges 
between -1 and 1. 



2.2. Numerical results on Joule heating 

Non- uniform heating and/or cooling profiles generate temperature gra- 
dients, even in good conductors such as copper. The magnitude of the effect 
is determined by the ratio of the energy flow due to electronic diffusion (ther- 
mal conductance) to that of the energy flow to the phonons. If diffusive flow 
is much larger, small temperature gradients exist and vice versa. Therefore, 
it is intuitively clear that more resistive samples have larger gradients, if 
the electron-phonon interaction does not depend on the mean free path I. 
However, as we discussed in the previous section, in impure metal films the 
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strength of the e-p interaction does actually depend on I, and a priori it is 
not fully clear how this affects the temperature profiles. 

To study the temperature profiles in mesoscopic metallic samples, we 
need to solve numerically the non-linear differential equation describing 
the heat flow in the system. We restrict ourselves to the simplest one- 
dimensional problem, which is valid for wires of approximately constant 
cross section. This is a good approximation for the samples described in 
the next section (Table H]) . The heat equation for a sample with uniform 
resistivity p reads 

Tx (j T(x) ^) = s[r(xr " T ^ ~ * t(x) + Ux) ' (3) 

where we have used the Wiedemann-Franz law for electronic thermal con- 
ductivity k = CT/p with C = Tr 2 kg/(3e 2 ) the Lorentz numberpland where 
qh(x) and q c {x) are the power density profiles for heating and cooling, re- 
spectively. If one heats the wire with a dc current density i, the Joule heating 
power density will be given by <jh(x) = pi 2 f(x), where f(x) = 1 at points 
where current flows, and f(x) = elsewhere. Typically, this problem has 
been solved for Dirichlet boundary conditions, for which temperature is fixed 
at the boundary. The Dirichlet problem describes well the case in which the 
wire is in direct contact with thick and wide normal electrodes,^ whose 
temperature stays constant. In contrast, our sample geometry is such that 
the wire does not have any contacts at the physical ends, and the heating 
current is passed through superconducting leads in direct contact with the 
wire (NS boundaries), so that no heat flow takes place through them. In 
this case, the correct boundary conditions are the von Neumann type, where 
dT/dx = 0. 

First, we discuss the numerical results on Joule heated wires. In this 
case the heated portion of the wire is much longer than the typical electron- 
electron scattering length at sub-Kelvin temperatures (L ~ 0.5 mm >> 
L e _ e ~ 1 pm), and therefore quasiequilibrium with a well defined electron 
temperature exists everywhere in the wire. In addition to the heated part of 
the wire, a short stub of length d = 9 pm (corresponding to Cu sample 2, see 
Table [1]) extends beyond the heated portion, where electrons can only diffuse 
and be cooled by phonons. This means that the temperature will start to 
fall, and the magnitude of this drop is a strong function of temperature itself. 
The first (and wrong, as we will see) intuition is that if the stub is much 
shorter than the length scale for phononic energy relaxation L e _ p = ^Dr e ^ p , 
where D is the diffusion constant, there is no significant effect below 1 K, as 
L e _ p ranges between ~1 mm and 10 pm at 100 mK -IK for typical thin 
films. 
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Fig. 1. Color online. Calculated temperature profiles with varying Joule 
power levels applied uniformly at positions x < 0. Solid (black) line uses 
material parameters for Cu sample 2 and dashed (red) line for the AIMn 
sample, see Tables [Q and [2] T p = 60 mK. 



Figure [T] shows the calculated temperature profiles for different Joule 
powers using the materials parameters for our experimental sample Cu 2 
(solid) and the AIMn sample (dashed) (Tables [T] and [5]) . The phonon tem- 
perature T p was set to 60 mK, which is the refrigerator temperature used in 
the experiments. The length of the stub was kept constant, using the value 
for the Cu 2 sample d = 9 p, The x-coordinate is such that Joule power 
is applied at x < 0, and the Joule power levels between Cu and AIMn are 
adjusted so that the electron temperatures in the bulk of the wires are equal. 
It is very clear that a short stub of this length does, in fact, have a signifi- 
cant effect on the profiles, contrary to our initial expectation. At T > 300 
mK, the temperature drop seems to be measurable for both materials, and 
stronger for AIMn, which has approximately five times higher resistivity, 
but also a weaker e-p scattering rate. Interestingly, the temperature drop 
extends mostly into the area x < 0, where Joule heat is being uniformly 
applied. Therefore, the bulk electron temperature, determined solely by the 



J.T. Karvonen, L.J. Taskinen and I.J. Maasilta 




-500 -400 -300 -200 -100 -500 -400 -300 -200 -100 



(a) Position (urn) (b) Position (|jm) 

Fig. 2. Calculated relative temperature profiles T/T^fc with varying Joule 
power levels applied at positions x < 0. Parameters used are the same as in 
Fig. [U (a) Cu sample 2, and (b) AIMn sample. Note that the curves have 
been shifted for clarity. Top curve corresponds to the lowest temperature 
Tbulk ~ 100 mK. 

electron-phonon scattering, can only be measured at > 40 /im distance away 
from the end of the wire. In the next section we describe the experiments 
used to study these temperature gradients. 

In Fig. [2j we also plot the same profiles scaled to the bulk values for 
each Joule power, so that the full temperature profile is clearly seen for all 
temperatures. We see that the length scale for the temperature drop grows 
strongly to > 100 /xm as one lowers the bulk temperature to 100 mK (top 
curve), although the relative drop becomes small. Also, the temperature 
profile is highly non-symmetric with respect to the average temperature, 
as expected for a non-linear system. At the higher temperature range, the 
AIMn profiles become clearly steeper due to the differences in p, £ and n. 

From the calculated profiles, one can see how the length scale and mag- 
nitude of the temperature drop depend on the bulk wire temperature. This 
is shown in Fig. [3j We have defined the energy loss length as the distance 
where T changes 90 % of the total change measured from the end of the stub. 
By comparing with the theoretical electron-phonon length L e ^ p in Fig. 3 (a), 
we see that our definition corresponds to roughly ~ 2.2L e _ p , with the cor- 
rect temperature dependence determined by the e-p scattering L e _ p ~ ynv 2 5 
where m = 3 for Cu and m = 4 for AIMn (Table [2|). The deviations at low 
temperatures are due to saturation caused by 2L e _ p approaching the length 
of the wire. The high temperature deviation in the case of AIMn shows that, 
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Fig. 3. (a) Symbols: Calculated energy loss length, defined as the length 
where T has changed 90 % of the total change. Squares: Cu, Triangles: 
AlMn. Solid lines show theoretical values 2.2L e _ p . Dashed lines shows the 
length of the wire, (b) Magnitude of the temperature drop, Squares: Cu, 
Triangles: AlMn. Solid lines are fits to power laws. Parameters used are the 
same as in Fig. [TJ 

at higher Joule power levels, the energy loss length is not simply given by 
~ 2L e _ p . We do not have a theoretical description for the magnitude of the 
temperature drop, which must be a function of the stub length d. However, 
we have simply plotted the values in Fig. [3] (b) with the observation that 
they follow a power law of the form AT ~ j i (Tt+2)/2 £ or materials, 
where (m + 2) and 2 are the exponents of the temperature dependencies of 
the heat flow rates due to e-p interaction and diffusion, respectively. 



2.3. Numerical results on NIS tunnel junction cooling 

So far we considered the case in which we had a non-uniform Joule 
heating profile. Technologically important is also the case in which a non- 
uniform cooling profile exists. This is in practice the case always for normal 
metal-insulator superconductor (NIS) tunnel junction electron coolers,^ as 
they never cover the whole area of the normal metal being cooled (there 
has to be room for thermometers). However, if the total dimensions of the 
normal metal island are small compared to L e _ p (regardless of how much 
area is uncooled), no significant gradients will develop. On the other hand, 
especially difficult is the case of phonon membrane coolers 041151 where the 
cooler junctions are located on the bulk of the wafer and extend a normal 
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Fig. 4. Color online. Calculated temperature profiles with a cooling power 
applied at positions x > 0. Solid (black) is the result for Cu sample 2 and 
(red) dashed line for the AIMn sample. Parameters used are the same as in 

Fig. m 

metal cold finger onto a thin insulating membrane. In that case, temperature 
gradients have to be taken into account. 

In Fig. 2] we plot the calculated temperature profiles for a long wire 
with a uniform cooling power applied at positions x > and for a constant 
phonon temperature T p = 340 mK, with Cu sample 2 and AIMn sample 
parameters (Tables[T][2]). The cooling power was chosen so that the minimum 
temperature would reach ~ 100 mK, which is a typical value for aluminum 
based tunnel junction coolers. ™ It is clear that the temperature starts to 
rise already within the cooled area, and rises back to the bulk value of T p 
within approximately 100 fj^rn, given by ~ 2L e _ p as before. This shows that 
effective cooling only works some tens of microns away from the junctions 
on the bulk substrate. When modeling membrane coolers more accurately, 
it is necessary to take into account the gradients of the phonon temperature 
^ and the fact that the electron-phonon interaction strength is weakened on 
membranes. 
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Fig. 5. (a) Schematic of the sample geometry and the measurement circuit. 
Black lines are the normal metal (Cu or AIMn), dark gray Al, and light gray 
Al or Nb. All the Al leads above the normal metal wire form NIS tunnel 
junctions, whereas the Al or Nb leads below have direct NS contacts. L is 
the length of the heated part of the wire and d is the length of the unheated 
part. d\ and d2 are the positions of the NIS-j unctions from the nearest NS- 
interface. (b) An SEM picture of an SN interface and a SINIS junction at 
the end of the heated wire of a representative sample. 



3. EXPERIMENTAL TECHNIQUES AND SAMPLES 

We performed non-uniform heating experiments on several Cu and AIMn 
wires. All samples were fabricated on oxidized or nitridized silicon chips 
by standard e-beam lithography and multi-angle shadow mask evaporation 
techniques. Figure [5] shows the schematic picture of the samples and the 
measurement circuit. Table [1] presents the essential dimensions of the sam- 
ples measured by SEM and AFM. The resistivity p was determined from the 
I-V measurement of the wire at 60 mK, from which the mean free path / was 
calculated using the Drude formula. 

All samples have a normal metal wire of length ~ 500 fim and width 
~ 400-600 nm, onto which two Al (Nb for Cu sample 2) leads form direct 
normal metal-superconductor (NS) interfaces. In addition, two pairs of Cu- 
AlO^-Al (NIS) tunnel junctions connect to the wire. The NS contacts are 
used to pass the heating current, and the SINIS junctions serve as electron 
thermometers in the middle and at the end of the wire. Near the heated 
wire, there is also a short, electrically isolated normal metal wire with a 
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SINIS thermometer, which measures the local phonon temperature. 

Table 1. Parameters of measured Cu and AIMn samples. L, d, d\ and efo 
are defined in Fig. [5j t is the thickness and A the cross sectional area of the 
normal metal wire, p is the measured resistivity, and I the mean free path. 
Parameter Cu sample 1 Cu sample 2 Cu sample 3 AIMn sample 



L [//m] 


473 


473 


492 


466 


d [pm\ 


11 


9 


20 


16 


d 1 [/xm] 


3 


2 


7 


7 


d 2 [/tm] 


8 


8 


14 


10 


t [nm] 


48 


32 


28 


55 


A [xlO" 14 m 2 ] 


1.5 


1.5 


2.5 


1.4 


p [xlO _8 Qm] 


2.5 


3.0 


3.8 


12.3 


I [nm] 


27 


22 


17 


3.2 



Because of Andreev reflection at the NS-j unctions, which are biased 
within the superconducting gap A, the Joule heating current in the normal 
metal is converted into a supercurrent in the superconductor, which does 
not carry any heat with it. Thus, the NS contacts are very good electrical 
conductors and, in the ideal case, perfect thermal insulators. This way the 
Joule heat does not leak into the superconducting side and the NS contacts 
do not cause any thermal gradients. In other words, the Joule heat is uniform 
between the the NS contacts. While the Joule current is being applied, we 
can simultaneously measure the applied power P from the measurement of 
current and voltage in four probe configuration, the electron temperature in 
the middle of the wire, where no gradients exist, as well as at the end of 
the wire, where significantly lower temperatures are expected. More details 
of the typical thermometer biasing and calibration are described in Ref. 17. 
The effect of the thermometers on the temperature profile was estimated to 
be insignificant. 



4. EXPERIMENTAL RESULTS 

As the thermometer approximately in the middle of the wire (the middle 
thermometer) measures the temperature in a region without thermal gradi- 
ents (Sec. 2.2), we see that Eq. ([3]) will reduce to = E(T™ — T"), and 
the values of n and S can be determined from the data. The thermometer 
at the unheated end of the wire (the side thermometer), on the other hand, 
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Fig. 6. Color online. Data from Cu sample 1. (a) The temperatures of 
the thermometers vs heating power density in log-log scale. Black line: 
experimental data from the middle thermometer. Gray line: experimental 
data from the side thermometer. Gray (Cyan) circles: numerical results for 
the middle thermometer, and gray (pink) stars: numerical results for the 
side thermometer. Dashed line is a guide to the eye T oc (P/V) 1 / 5 . (b) The 
logarithmic numerical derivatives of the experimental middle thermometer 
data, and (c) the same for the side thermometer. 



lies in the region of a gradient. As it comprises of two NIS-tunnel junctions 
separated by a small, but a significant distance, the two junctions measure 
two different temperatures, and the combined SINIS measurement will give 
a value in between these two temperatures. The measured temperature is 
not necessarily the average of the two due to the nonlinearity of the NIS 
thermometer. 

Cu samples 1 and 2 have a similar sample geometry, except for a differ- 
ence in the thickness t, and in the electron mean free path I. In Fig.s[6ja) 
and [7|[a) we plot the temperatures of both thermometers vs heating power 
density P/V in log-log scale. Clear temperature difference between the two 
thermometers can be seen at P/V > 1 pW/fim 3 , and at 100 pW///m 3 the 
difference is roughly 50 mK. Below 0.1 pW//im 3 both temperatures saturate 
mostly because of noise heating. The theoretical points have been calcu- 
lated by solving Eq. ([3]) using the appropriate sample parameters. In the 
calculation, we used values for T, obtained from the middle thermometer 
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Fig. 7. Data from Cu sample 2. Symbols are explained in the caption of El 



data (Table EJ). The temperature of the side SINIS thermometer is defined 
in our analysis as the average of the two NIS-junction temperatures. This 
definition is justified, because the difference between the two temperatures 
is small and within the size of the plotted datapoints for all sample geome- 
tries studied in this paper. From Figs.s Eta) and [T^a) it is clear that the 
theory agrees well with the experimental data, showing that the observed 
temperature difference is fully explained by phonon cooling and diffusion. 

From the phonon thermometer data (not shown), we have seen that 
T e 3> T p for all samples in this work, and therefore for the middle ther- 
mometer we can approximate Pheat/V = ST™, where V is the volume of the 
heated portion of the wire. The temperature dependence n and strength E 
of electron-phonon interaction can thus be obtained from the plots of the 
middle thermometer temperature versus heating power density in s[6ta) and 
[7|(a). A more detailed look at n can be obtained by plotting the logarithmic 
derivatives d(logT e )/d(log P) = 1/n, shown in Figs.s EKb) and[7Jb). Mea- 
sured data in both samples scales clearly as P oc T e 5 . However, because of the 
temperature drop at the end of the wire, the data from the side thermome- 
ter, Figs. [61(c) and E|c), show a temperature dependence P oc T e 5,5 . This 
exponent does not correspond to the actual power law of the e-p interaction. 

Copper sample 3 is a bit thinner and has a longer unheated end section 
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Fig. 8. Data from Cu sample 3. Symbols are explained in the caption of 
Fig. B 



compared to samples 1 and 2 (Table [T]). From Fig. Eta) we observe that 
the measured temperature difference between the two thermometers is much 
larger than for samples 1 and 2, and also larger than what the theoretical 
calculation predicts. We do not fully understand this at the moment, but it 
is possible that the Lorentz number is reduced due to inelastic scattering on 
the surface.^ Indeed, high resolution SEM images of sample 3 showed that 
the surface was more irregular. Nevertheless, we can obtain the temperature 
dependence of e-p interaction from the middle thermometer [Figs. Eta), (b)], 
showing again an agreement with P oc Tj* , while the data from the side 
thermometer fits P oc Tj. 

The last sample is made from aluminum doped with 0.7 at % manganese. 
Due to the high impurity concentration, the mean free path I is much shorter 
than for copper wires, and the AIMn wire is ~ 4 times more resistive (Table 
[I]). Figure Eta) shows that the measured data and the numerical result agree 
well. It may be surprising, perhaps, that the temperature difference between 
the thermometers is actually smaller than in Cu sample 3. Although the 
AIMn sample is more resistive and thus diffusion is weaker, the e-p scattering 
is much weaker. Therefore, the unheated end of the wire is not as effectively 
cooled by phonons in AIMn. 

The sample is clearly in the limit < 1 and the scattering potential is 
dominated by the Mn impurities. The data from the middle thermometer, 
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Fig. 9. Data from the AIMn sample. Symbols are explained in the caption 
of Fig. El 



Figs.s[9^a) and (b) show that P oc Tg . This is in agreement with the theory 
including vibrating scatterers, ^ and has also been observed for other Mn 
concentrations.^ Again, the side thermometer gives a much higher apparent 
temperature dependence: P oc Tj. 

The values for coupling constants £ can be determined from the slopes 
of the (P/V, T)-plots in logarithmic scale as fitting parameters, keeping n 
fixed at n = 5 or n = 6. In Table [5] the measured values for £ are pre- 
sented for all the samples. For Cu samples 1-3, S decreases as a function 
of electron mean free path I; in other words, the electron-phonon interac- 
tion weakens with increased purity of the samples. This is evidence that 
the theory for pure electron-phonon coupling does not apply in our Cu thin 
films, although the temperature dependence agrees with the simplest theory 
without disorder Possible explanations for the experimental result P oc T 5 
are that our Cu samples are either in the transition region ql ~ 1 or that 
the scattering potential is not fully vibrating. Our earlier conclusions^^ on 
the temperature dependence of the electron-phonon coupling in disordered 
Cu and Au films were not correct, because the temperatures were measured 
with a side thermometer. For aluminum manganese samples with varying 
impurity concentration, S is linearly dependent on I, consistent with the 
fully vibrating disorder theory!^ 
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Table 2. Measured values for coupling constant S 

Cu sample 1 Cu sample 2 Cu sample 3 AIMn sample 

n 5 5 5 6 

E [W/K"m 3 ] 1.8 xlO 9 2.1 xlO 9 2.5 xlO 9 3.4 xlO 9 

5. CONCLUSIONS 

We have shown that thermal gradients are easily generated in meso- 
scopic samples at sub-Kelvin temperatures, even for good conductors such 
as copper. This fact has a strong effect on studies of thermal properties and 
thermometry. To obtain correct information on the electron-phonon interac- 
tion strength, for example, one has to make sure that electron temperature 
is measured in a location without thermal gradients. In addition, our re- 
sults also imply that for tunnel junction coolers, large cold fingers outside 
the junction area are not effectively cooled. Also, even if the non-heated or 
non-cooled area is smaller than the electron-phonon scattering length -L e _ p , 
thermal gradients will develop as long as the total size of the normal metal 
is of the order of L e _ p . 
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